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Abstract 

In this paper, stabihty theorems for stochastic differential equations and back- 
ward stochastic differential equations driven by G-Brownian motion are obtained. We 
show the existence and uniqueness of solutions to forward-backward stochastic differen- 
tial equations driven by G-Brownian motion. Stability theorem for forward-backward 
stochastic differential equations driven by G-Brownian motion is also presented. 
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1 Introduction 

Consider a family of ordinary stochastic differential equations (SDEs for short) parameterized by 
e > 0, 

Xf = xl + [ b'{s,XI)ds+ [ a'{s,XI)dWs,te[0,T], 
Jo Jo 

where Wt is classical Brownian motion. It is well known that the strong convergence of the coefficient 
in implies the strong convergence of the solutions, that is, if 

and 

E[ffi\¥is,X^) - + \a%s,X^) - a%s,X^,)\^)ds] ^ 0, as 6 ^ 0, 

then under Lipschitz and other reasonable assumptions, their solutions also converge strongly in L^, 

V t G [0, T],E[\Xf - X^op] ^ 0, as e ^ 0. 
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This result, known as the continuous dependence theorem, or the stabihty property, can be found 
in many standard textbooks of SDEs (e.g., see [16]). 

Backward stochastic differential equations (BSDEs for short) driven by classical Brownian motion 
were introduced, in linear case, by Bismut [3j in 1973. In 1990, Pardoux and Peng considered general 
BSDEs (see [12]). Similar continuous dependence theorem for the case of backward stochastic 
differential equations was obtained by El Karoui, Peng and Quenez (1994) [6] and Hu and Peng 
(1997) do]. 

As for the forward-backward equations, Antonelli [Ij first studied these equations, and he gave 
the existence and uniqueness when the time duration T is sufficiently small. Using a PDE approach. 
Ma, Protter and Yong [llj gave the existence and uniqueness to a class of forward-backward SDEs 
in which the forward SDE is non-degenerate. In 1995, Hu and Peng [9j study the existence and 
uniqueness of the solutions to forward-backward stochastic differential equations without the non- 
degeneracy condition. 

Motivated by uncertainty problems, risk measures and the superhedging in finance, Peng(2006, 
see |13] ) has introduced the notion of sublinear expectation space, which is a generalization of 
classical probability space. Together with the notion of sublinear expectation, Peng also introduced 
the related G-normal distribution and G-Brownian motion. The expectation associated with G- 
Brownian motion is a sublinear expectation which is called G-expectation. The stochastic calculus 
with respect to the G-Brownian motion has been established by Peng in [13] . |14j and [15]. Since 
these notions were introduced, many properties of G-Brownian motion have been studied by authors, 
for example, [5], [7], [8], [I7]-[l9], et al. 

Therefore, the natural questions are: stability properties for stochastic differential equations 
and backward stochastic differential equations driven by G-Brownian motion are also true? How 
to obtain the existence and uniqueness of the solution of a forward-backward stochastic differential 
equations driven by G-Brownian motion? The goal of this paper is to study stability properties for 
stochastic differential equations driven by G-Brownian motion (G-SDEs for short) and backward 
stochastic differential equations driven by G-Brownian motion (G-BSDEs for short). Indeed, under 
Lipschitz or integral-Lipschitz condition and other reasonable assumptions, stability theorems for 
G-SDEs and G-BSDEs are obtained. Meanwhile, we also show the existence and uniqueness of the 
solution of a new type of forward-backward stochastic differential equations driven by G-Brownian 
motion. 

This paper is organized as follows: in Section 2, we recall briefly some notions and properties 
about G-expectation and G-Brownian motion. In Section 3, we study the stability properties of 
G-SDEs, while Section 4, study the G-BSDEs case. At last, the existence and uniqueness of the 
solution of forward-backward stochastic differential equations driven by G-Brownian motion are 
obtained. Stability theorem for forward-backward stochastic differential equations driven by G- 
Brownian motion is also presented. 

2 Preliminaries 

In this section, we introduce some notations and preliminaries about sublinear expectations and 
G-Brownian motion, which will be needed in what follows. More details concerning this section may 
be found in [l3], [H] and [E]. 

Let be a given set and let 7i he a, linear space of real valued bounded functions defined on Q. 
We suppose that 71 satisfies C £ Ti ior each constant C and \X\ £ T-l, if X £ 7i. 
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Definition 2.1. A sublinear expectation E is a functional E : H — >■ -R satisfying 

(i) Monotonicity: E[X] > E[Y] iiX>Y. 

(ii) Constant preserving: E[C] = C for C G R. 

(iii) Sub-additivity: For each X,Y eH, E[X + ¥]< E[X] + E[Y]. 

(iv) Positive homogeneity: E[XX] = XE[X] for A > 0. 

The triple {Q,T-L,E) is called a sublinear expectation space. If (i) and (ii) are satisfied, E[-] is 
called a nonlinear expectation and the triple (r2,'H,E) is called a nonlinear expectation space. 

Prom now on, we consider the following sublinear expectation space (r2,^,E): ifXi,--- ,X„ eV., 
then (p{Xi,- ■ ■ , Xn) G for each (p G Ci^upO^"'), where Ci^Lip{R^) denotes the linear space of func- 
tions ip satisfying \p{x) — ^{y)\ < C(l + \xY^ + |y|™')|a; — y\ for x,y e i?", some C > 0,m G N 
depending on (p. 

Definition 2.2 . Let X and Y be two n-dimensional random vectors defined on nonlinear expec- 
tation spaces (f2i,Hi,Ei) and {^2,^2,^2), respectively. They are called identically distributed, 
denoted hy X = Y, if 

Ei[v?(X)] = E2[ip{Y)], for V99 G Q,Lip(i?"). 

Definition 2.3. In a nonlinear expectation space (f2,'H,E), a random vector Y € "H" is said to be 
independent from another random vector X G T-L"^ under E[-], if 

E[p{X,Y)] = E[E[^{x,Y)]^=x], for G 

X is called an independent copy oi X ii X = X and X is independent from X. 

Definition 2.4 (G-normal distribution). In a sublinear expectation space (0,H,E), a random 

variable X eH with 

E[X^] = a^,-E[-X^] =a^, 
is said to be A^(0; [ct^, CT^])-distributed, if for each X eH which is an independent copy of X we have 

aX + bX = + b'^X, Va, b>0. 

Definition 2.5 (G-Brownian motion). A process {Bt{u!)}t>o in a sublinear expectation space 
(n, H, E), is called a G-Brownian motion if for each n G iV and < ii < • • • < t„ < 00, Bt^ , • " " > -St„ G 
^ and the following properties are satisfied: 

(i) Bo(w) = 0; 

(ii) Por each t,s>0, the increment -Bj+s — Bf is A^(0; [a^, (7^])-distributed and is independent from 

{Bt^ ,■■ ,BtJ for each nE N and < ti < ■ ■ ■ < < t. 

We denote by $7 = Cq{R'^) the space of all i^'^-valued continuous paths (wt)tgfl+, with ojq = 0, 

00 

equipped with the distance p{oj^,uj'^) := J2 2~*[(max \ujj — |) A 1]. Considering the canonical 

i=i t&[o,i] 

process Bt{uj) = {uJt)t>o- For each fixed T > 0, set := {w.at oj e CI} and 

LipiClT) ■■= MBt„Bt„...,BtJ : m > G [0,r],^ G Q,Lip(i?'^^™)}, 

and define Lip(Q) := (J Ljp(fi„). 

n=l 
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Let ^ be a G-normal distributed, or A7^(0; [cr^, l])-distributed random variable in a sublinear 
expectation space (r2,^,E). We now introduce a sublinear expectation E defined on Lip(yi) via the 
following procedure: for each X G Lip{Q) with 

X = (p{Bt^ - Bto,Bt2 - Bti,--- , Btrr, - -St^-l)) 
for some cp G Ci,Lip(-R'^^"*) and = to <*!<•••< < OO) we set 

E[v{Bt, -Bt„--- ,Bt^- Bt^_,)] := E[<p{Vh^^i, ■■■ , V^m - im-iCm)], 

where (^i, • • • ,S,m) is an m-dimensional G-normal distributed random vector in a sublinear expecta- 
tion space (Q,,7i,E) such that = A^(0; [cr^, 1]) and such that ^j+i is independent from (,^1, • • • 
for each i = 1, - ■ ■ ,m. 

The related conditional expectation of X = ip{Bt^ — Bt^^Bt^ — Bt^, ■■■,Bt^ — Bt^_^) under VLt- 
is defined by 

t[X\nt^]=n^{Bt,,B,,-Bt,,... ,Bt^-Bt^_,)\^t^] 
■.= iP{Bt„,--- ,Bt.-Bt^_^), 

where 

tp{xi, ■■■ ,Xj)= E[(p{xi, ■ ■ ■ ,Xj, ^ytj+l - tj^j+i, ■■■ , ^Jtm- tm-\im)\- 

Definition 2.6. The expectation ]E[-] : Ljp(O) R defined through the above procedure is called 
G-cxpcctation. The corresponding canonical process {Bt)t>o in the sublinear expectation space 

(O, Ljp(r2), E) is called a G-Brownian motion. 

We denote by Lq(^}t),p > 1, the completion of Lip^nx) under the norm \\X\\p := (]E[|X|^])-^/^'. 
Similarly, denote Lq{^) is complete space of Lip{U). We give some important properties about 
conditional G-expectation E[-|r2t],t G [0,r]. 

Proposition 2.1. The conditional expectation E[-|r2t],t G [0,r] holds for each X,Y E L(j(f2t) : 

(i) If X > y, then E[X|Ot] > E[y|Jlt]. 

(ii) E[77|r2t] = 77, for each t G [0,oo) and ij G L^(Qt). 

(iii) E[X\nt] - E[Y\nt\ < E[X - Ypt]. 

(iv) E[r]X\nt] = r]+E[X\nt] + r]-E[-X\nt] for each bounded rj G ^^.(nt). 

(v) E[E[X|f]t]|Jls] = E[X\ntr,s], in particular, ]E[E[X|Jlt]] = E[X]. 

Next, we introduce the Ito's integral with G-Brownian motion. For T G R^, a partition ttt of 
[0, r] is a finite ordered subset ttt = {to,ti, ■■■,tN} such that = to < < •■• < ^at = T, 

^(ttt) ■■= max{|tj+i - ti\ : i = 0, 1, N - 1}. 

Using vr-^ = {t^, , t^} to denote a sequence of partitions of [0, T] such that lim fJ-{iT^) = 0. 

Let p > 1 be fixed. We consider the following type of simple processes: for a given partition 

N-l 

TTT = {to,ti,...,tN} of [0,T], set r]t{uj) = Yl 6('^)I[tfc,tfc+i)(*), where ^a; e L^{^tk),k = 0, l,...,iV-l 

k=0 

are given. The collection of these processes is denoted by Mq°(0,T). For each p > 1, we denote by 
Mg([0,r];i?") the completion of Mg'°([0, T]; i?") under the norm \\vt\\M^{[o,T]) ■= {Jo ^ihl^ldi)^ ■ 
Definition 2.7. For an 77 G Mq°(0,T), the related Bochner integral is 

T N-1 

/ r]t{u)dt := (,k{i^){tk+i - tk). 
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Let {Bt)t>o be a 1-dimensional G-Brownian motion with G{a) := ^W\_aB\\ = ^(c^a^ — g^a ), 
where = E[Bl],g^ = -t[-Bl], < a < a < oo. 

Definition 2.8. For an r/ e M^'"(0,r) of the form rit{uj) = E Cfc(^^)I[tfc,tfc+i)(*)> define 



A;=0 

iV-1 



/ r^{s)dB,:=^^k{Bt,^,-Bt,). 



Proposition 2.2. For each r] G Mg'°(0,r), then 

E[/ r?(s)(iS,] = 0, E[(/ r]is)dBsf]<a^ [ E[r]\s)]ds. 
Jo Jo Jo 

Definition 2.9. For the 1-dimensional G-Brownian motion Bt, we denote {B)t is the quadratic 



variation process of Bt, where {B)t := lim J2 (-^t^ ~ B^n)'^ = (Bt)^ ~ Jo BgdBg. 

/i(7rf )^0 k=0 ''+^ 

Definition 2.10. For each r/ G M^'°(0,T), define r]{s)d{B)s := E ^fe((^)tfc+i - (^)tj- 

fe=o 

Proposition 2.3. For any < t < T < oo, 

(i) E[| ffvtd{B)t\] < a^E[J^ \rjt\dt],y vt G M^{0,T). 

(ii) ]E[(/o^77idi?i)2] = E[/„^r,2d(B),],V Tit G M2(0,r). 

(iii) E[/(f \rjt\Pdt] < J^^ E[\rjt\P]dt,y r^t G Mg(0,r), p > 1. 

3 Stability theorem of G-stochastic differential equa- 
tions 

In this section, we consider the stability theorem of G-stochastic differential equations. Consider the 
following stochastic differential equations driven by d-dimensional G-Brownian motion: 

ft ft ft 

Xt = Xo+ b{s,Xs)ds+y] / hij{s,Xs)d{B\B^), + Y. X,)dBi, t e [0,T], (3.1) 

Jo ij^^Jo ^^-^Jo 

the initial condition Xq G i?", and b,hij,aj are given functions satisfying b{-,x),hij{-,x), aj{-,x) G 
M^([0, T];R^) for each x G i?". Consider the following G-SDEs depending on a parameter £(e > 0): 

ft ft ft 

XI = X^+ b%s,Xl)ds+Y^ / Kj{s,Xl)d{B\B^)s + Y, cj'j{s,Xl)dBl iG[0,r]. (3.2) 

We make the following assumptions: 
Assumption 3.1. For any e > 0,x G i?", 6^(-, x), /ifj(-, x), cr|(-, x) G M^([0, T]; i?'*), G i?". 
Assumption 3.2. For any e > 0, x, xi, a;2 G -R" : 

(HI) |6^(i,x)p+ E |/i|,(i,x)p+Ek|(i,x)p<a2(i) + «2(i)|^|2^ 
i,j=i j=i 
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(H2) \b%t,xi)-b'{t,X2)\^+ E \hl^{t,xi)-hl.{t,X2)\HY: k^^(t,xi)-a|(t,X2)|2 < a^{t)p{\xi-X2\^), 

i,j=l j=l 

where ai G M^{[0,T]),a2 ■ [0,T] — )• and a : [0,T] — i2+ are Lebesgue integrable, and p : 
(0,+oo) — )• (0, +00) is continuous, increasing, concave function satisfying p(0+) = 0, Jq '^dr = 

+CO. 

Assumption 3.3. (i) V t G [0,r], as e ^ 0, 

where (p = b, hij and aj, respectively, i, j = 1, • • • ,d. 
(ii) As e 0, 

XI ^ AO. 

Remark 3.1. The Assumptions 3.1 and 3.2 guarantee, for any £ > 0, the existence of a unique 
solution Af e M2([0,r];i?") of G-SDEs (3.2)(see pj), while the Assumption 3.3 will allow us to 
deduce the following stability theorem for G-SDEs. 

Theorem 3.1. Under the Assumptions 3.1, 3.2 and 3.3, we have the following convergence: as 
e ^ 0, 

V t G [0,r], E[\Xf - Af |2] ^ 0. (3.3) 

In order to prove Theorem 3.1, we need the following lemmas: 
Lemma 3.1 (see Chemin and Lerner [4]). Let p : (0, +00) — )• (0,+oo) be a continuous, increasing 

h p{r)' 

defined on (0, +00) satisfying 



function satisfying p(0+) = 0, 'TP)df = +c« and let u be a measurable, nonnegative function 



u{t)<a+ I P{s)p{u{s))ds, tG(0, +00), 
^0 



where a G [0, +00), and /3 : [0,T] — )■ is Lebesgue integrable. Then 

(i) if a = 0, then u{t) = 0, for t £ [0, +00); 

(ii) if a > 0, then 

uit) < u~\t;(a) + / I3{s)ds), 



where v{t) := //^ ■^ds,to G (0,+oo). 

Lemma 3.2 (see Peng [151). Let p : R ^ R he a continuous increasing, concave function defined 
on R, then for each A G L^(Q),V t > 0, the following Jensen inequality holds: 

p{E[x\nt])>K[p{x)\nt]. 

Proof of Theorem 3.1. Let Af := Af - A^, Ag := Ag - A^, then 



Af = Ag + (\b^{s,Xl) - b\s,X^,))ds 
Jo 

+ E l\KMXl)-h%{s,X^^))d{B\B^) 

+ E [\<^j{s,Xl)-a^{s,X''j)dBi, 
j=i 



(3.4) 
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and 

\Xt\^ < + I [\b'{s,Xl) - b%s,X^,))ds\' + I j\h'{s,Xt) - b\s,X^,))ds\^ 

Jo Jo 

+ E I [\hUs,XI)-ht^{s,Xf))d{B\B^)s\'+Yl I [\hlj{s,X^^)-h%{s,X'j)d{B\B^)s\' 
+ El Aa|(.,Xf)-a|(.,xO))di?i|2 + ^| Aa|(.,xO)-aO(.,xO))d5i|2}, 

(3.5) 

taking the G-expectation on both sides of the above relation and from Proposition 2.3, we get 

n\X^f] < C{|lgp + fn\b%s,Xl) - b%s,X^)f]ds + fE[\b%s,X',) - b\s,X',)f]ds 
Jo Jo 

+ E j\\hUs,Xl) - hl^{s,X^^)\^]ds + j\\hl^{s,X^^) - h%{s,X^^)\^]ds (3_g) 

+ fE[\<^'jis,Xl)-a'^{s,xX]ds + Yl f E[\a^s,X',) - a^{s,X^,)\']ds}, 

by Assumption 3.2, we have 

E[\X^\'] < C%T) + C2 f a\s)E[p{\XI\')]ds, (3.7) 

Jo 



where 



C%t) : = C fEWis,X^,) - b'{s,X',)\ys + CY^ f nKjis^Xt) - h%{s , X^,)\^\ds 
Jo ,^^^,Jo 

+ CJ2 fnW'j{s,X'j - a^{s,X^^)\']ds + C\XI\'. 



Because p is concave and increasing, from Lemma 3.2, we have 

E[\Xn''] < C%T) + C2 f a\s)piE[\Xlf])ds. (3.8) 

Jo 

Since as e — >■ 0, C^{T) 0, hence, from Lemma 3.1, we get 

]E[|X||2] ^ 0, as e ^ 0. 

The proof is complete. 

A special case of Assumption 3.2 is 
Assumption 3.4 (Lipschitz condition). For any x\,X2 G BP', there exist constant Co > such 
that 

|#(t,,xi) -(/.^(t,X2)| <Co|xi-X2|, [0,r], 
where 4> = b, hij and aj, respectively, i,j = l,--- ,d. 

Corollary 3.1. Under the Assumptions 3.1, 3.3 and 3.4, we have the convergence of the solution of 
the G-SDEs (3.2) in the sense of (3.3). 
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4 Stability theorem of G-backward stochastic differen- 
tial equations 

In this section, we give a stability theorem of backward stochastic differential equations driven by 
d-dimensional G-Brownian motion (G-BSDEs for short). Consider the following type of G-backward 
stochastic differential equations depending on a parameter (5 > 0) : 

Yf = E[C'+ r f{s,Y^)ds+ V r gf^{s,Y,')d{B\B^)s\nt], t e [0,T], (4.1) 

i,j=i ''^ 

where £ Ll{aT;R^) is given, and f\-,y),gf^{;y) £ M^{0,T;R^). 

We further make the following assumptions: 
Assumption 4.1. For any 6 > 0,y,yi,y2 G -R", 

(HI) \f{t,y)\+ E |4.(t,y)| </3(t) + C|y|, 

m\f{t,y,)-f{t,y2)\+ E \g!,{t,y,)-gf^{t,y2)\<p{\yi-y2\), 

where C > 0,/3 G M^([0, T]; i?+), and p : (0, +oo) — )■ (0, +oo) is continuous, increasing, concave 
function satisfying p(0+) = 0, Jq -^dr = +oo. 
Assumption 4.2. (i) V i G [0,r], as 5 ^ 0, 



where </> = /, gij respectively, j = 1, • • • , d. 
(ii) As (5 ^ 0, 

E[ie'-c°|]^o. 

Remark 4.1. Under the Assumptions 4.1 and 4.2, G-BSDEs (4.1) has a unique solution. The proof 
goes in a similar way as that in [2], and we omit it. 

Theorem 4.1. Under the Assumptions 4.1 and 4.2, we have the following convergence: as 5 — t- 0, 

vte[o,r], E[|y/ - y^^i] ^ 0. (4.2) 

Proof. Let y/ := - Y^,^^ := - then 



\Y,'\<m'\+ [ \f'{s,Y,')-f{s,Y,^)\ds+Y. I \9Us,Y!)-g%{s,Y!)\d{B\B^),\nt] 

< E[|e^| + \f{sX) - f{s,Y!)\ds + f2 [ - g%{sXmB\B^), (4.3) 



d „y 



+ / \f\s,Y!)-f{sX)\ds+Y, I \gtM'^!)-gt,{s.Y^)\d{B\B^),m. 
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Taking the G-expectation on both sides of (j4.3p . we have 
E[|y/|] < m'\] + rmHs,Y,')-f{s,Y,^)\]ds + Cj2 r n\gU',Y,') - g^^{s,Y,')\]ds 

+ / t[\f{s,Y!)-f{s,Y!)\]ds + CY,f nij{s,Y^)-gt^{s,Y^)\]ds. 
From the Assumption 4.1, Propositions 2.1 and 2.3 as well as Lemma 3.2, we have 



nyt\\<c'{Q)+K^ np{\Y^\)]ds 
<c\q)+k^ j\my!\\)ds. 



(4.5) 



where 

C^(0):=E[|f|]+ r nf{sX)- f{s,Y!)\]ds + CY, r E[\9Us,Y,') - 9U',ys)\]ds. 
Jo ^j^^Jo 

Since as 5 — )■ 0, C^(0) — >• 0, hence, from Lemma 3.1, we have 

]E[|y/|] ^ 0. 

The proof is complete. 

A special case of Assumption 4.1 is 
Assumption 4.3. For any S > 0,yi,y2 & R", there exist constant Cq > such that 

\(t>'{t,yi)-4>Ht,y2)\<Co\yi-y2\, te[0,T], 

4> = f,9ij respectively, i,j = !,■■■ ,d. 

Corollary 4.1. Under the Assumptions 4.2 and 4.3, we have the convergence of the solution of the 
G-BSDEs (4.1) in the sense of (4.2). 



5 Forward- backward stochastic differential equations 

The goal of this section is to show the existence and uniqueness of forward-backward stochastic 
differential equations driven by G-Brownian motion. For notational simplification, we only consider 
the case of 1-dimensional G-Brownian motion. However, our method can be easily extend to the 
case of multi-dimensional G-Brownian motion. We consider the following system: 

Xt = x+ I bis,Xs,Ys)ds+ I h{s,Xs,Ys)d{B)s+ [ a{s,Xs,Y,)dB„ 

Jo Jo Jo i\ 

< rp rp (5.1) 

Yt=E[^+ [ fis,Xs,Ys)ds+ [ gis,Xs,Y,)d{B)s\nt], tG[0,T], 
Jt Jt 

where the initial condition x G i?, the terminal data ^ G Lq{^t'-, R), and b,h,a,f,g are given 
functions satisfying b{-,x,y),h{-,x,y), a{-,x,y), f{-,x,y),g{-,x,y) G M^{[0,T];R) for any {x,y) E 
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(5.2) 



and the Lipschitz condition, i.e., \(j){t,x,y) — (j){t,x' ,y')\ < K{\x — x'\ + \y — y'\), for each t E 
[0, T], (x, y) G R^, {x' , y') E R^, cj) = b,h,a, f and g, respectively. The solution is a pair of processes 
{X,Y) e M2(0,r;i?) X M^{0,T;R). 

This model is called forward-backward because the two components in the system ()5.ip are 
solutions, respectively, of a G-forward and a G-backward stochastic differential equation. 

We first introduce the following mappings on a fixed interval [0, T] : 

A* : M^(0,T;R) x M^{0,T;R) M^{0,T;R) x M^{0,T;R),i = 1,2, 
by setting Aj,i = l,2,t £ [0,T], with 

A]iX,Y) = x+ [ b{s,Xs,Ys)ds+ [ h{s,Xs,Ys)d{B),+ [ a{s, Xs,Ys)dBs, 
Jo Jo Jo 

A2(X,y) =E[^ + ^ f{s,Xs,Ys)ds + g{s,Xs,Ys)d{B)s\nt]. 

Lemma 5.1. For any (X, Y), {X',Y') G M^(0, T; R) x M^(0, T; i?), we have the following estimates: 

n\Aj{X,Y)-Aj{X',Y')\^]<C ['E[\Xs-X'f + \Ys-Y:\^]ds,te[0,T], 

Jo 

E[|A2(x,y) - A2(x',y')l'] < C j\[\Xs-x'f + \Ys-Y^\']ds,t G [o,r], 

where C = 24K'^,C' = SK"^, K is Lipschitz coefficient. 
Proof. ^ 

E[|AKx,y)-A,\x',y')l']<4 / E[|6(5,x„y,)-6(s,x^,yi)|2]ds 

JO 

+ 4 / E[|/i(s,x„y,)-/i(s,x;,y;)|2]ds 

JO 

+ 4 / E[|a(s,X„y,)-a(s,X;,y;)|2]ds 
Jo 

|A?(X,y) - A?(X',y')|2 < 2E[| £ fis,Xs,Ys) - fis,X',X)ds\^ 

+ I J\{s, Xs,Ys) - gis, X',X)ds\''\^t], 



(5.3) 



And since 
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then 



n\AUx, Y) - Af{x', y )p] < 2E[| J^^ f{s, Xs,Y,) - f{s, x^, y;)dsp 

+ \ j\is,Xs,Ys) - gis,XiX)ds\^] 

< 2 E[|/(s, Xs,Ys) - f{s, X'„ 
+ 2 EMs, Xs,Ys) - gis, X'^, 

< 8K^ J E[\Xs -XX + \Ys - 

Let us consider the space M^(0, T; R) x Mq{0, T; R), with the norm \ \{X, Im^{o,t)xM^(o,t) '•— 

ll^llAf2(o,T) + II^IIm2(o,t) = Io^[\Xs\'^]ds + J^E[\Ys\'^]ds, this is a Banach space. 
Theorem 5.1. Under the previously stated hypotheses and under the additional assumption that 
{2^/6 + 2V2)KVt < 1, then there exists a unique solution {X,Y) G M^{0,T;R) x M^{0,T;R) of 
the forward-backward stochastic differential equation (|5.ip in the M^(0,T) x Af^(0, T) sense. 
Proof. Let us consider the space M^{0,T;R) X M^{0,T;R), with the norm 

II(-'^>"^)IIm2(0,T)xM2{0,T) ■= II-'^IIm2{0,T) + II^IIm2(o,t)- 

/ X^(x Y) \ 

We can view the system (15. ID as the operator At{X,Y) := ( \^(^x y) ) ' ^^^^ 



||At(x,y)-Ai(x',y')llM2(o,T) xM2{0,T) 

= ||AKx,y) - Ai(x',y')llM2(o,T) + ||A2(x,y) - A2(x',y')llM2(o,T) 



= (/ E[(Ai(x,y)-AKx',y'))']^i)^ (^-^^ 

JO 

+ ( rE[(A?(x,y) - A?(X',y'))']di)^- 

JO 

Because the terminal data ^ G Lq{Qt', R), and 6, /i, a, f, g are given functions satisfying a;, y), x, 
a{-,x,y), f{-,x,y),g{-,x,y) G Mg.([0, T]; i?) for any (x,y) G i?^ and the Lipschitz condition, we can 
prove ||At(X,y)||j^2(o,T)xM2(o,r) < +oo,V(X,y) G M^{0,T;R) x M2(0,r;i?). Next, we prove it 
is a contraction mapping. From the Lemma 5.1, we can obtain 

\\At{X,Y) - At(X',y')||Af2(0,T)xM2(0,T) 
rT ft 



<{f C [ E[\Xs - x'f + |y, - Y:\^]dsdt)-2 

Jo Jo 

+ {[ C f E[|x, -x^|2 + |y,-y;|2]dsdt)5 (5.5) 

Jo Jt 

fT 

<{VC + VC^)Vt E[\Xs-X'J^ + \Ys-Y;\^]ds 
Jo 



{Vc + Vc^)VT\\ix -x',Y- y')||M2 



(0,r)xM2 (0,T)- 
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From the assumption (2\/6 + 2\/2)^'^\/r < 1, we can obtain that Aj(X, Y) is a contraction mapping. 
Hence a unique fixed point for A exists and this is the solution of our system (|5.ip . The proof is 
complete. 

In the last section, we present stability theorem for forward-backward stochastic differential 
equations driven by G-Brownian motion. 

6 Stability theorem of forward- backward stochastic dif- 
ferential equations 

Consider a family of forward-backward stochastic differential equations with parameter (7 > 0), 

X^ = x^+ fh"'{s,X],Yr)ds+ f h"'{s,X],Yr)d{B)s+ f a\s,X2,Y^)dBs, 

^° ^ ^ (6.1) 

Y;< = ne<+f P{s,X],Y^)ds+ f g\s,X],Y^)d{B)s\ntlte[f),Tl 
Jt Jt 

where the initial condition x'^ G R, the terminal data G Lq{Qt', R), and IP ,a'^ , are 
given functions satisfying b'^{-,x,y),h'^{-,x,y), a'^{-,x,y), p{-,x,y),g^{-,x,y) G M^{[0,T];R) for 
any (x,y) G R^ and the Lipschitz condition, i.e.. 
Assumption 6.1. 

\(t>{t, X, y) - (t>{t, x', y')\ < Ki\x - x'\ + \y-y'\), 

for each t G [0, T], (x, y) G R^, {x' , y') G R^, (j) = b'^ ■: h"^ 1 P and ^ respectively. 

We further make the following assumption: 
Assumption 6.2. (i) V t G [0,r], as 7 0, 

* tW{s, Xl - Xl Y,^)\^]ds ^ 0, 





where (j) = b,h,a, f and g, respectively, 
(ii) As 7 0, 

x^^x°, E[\C-f\^]^0. 

Theorem 6.1. Under the Assumptions 6.1 and 6.2, then as 7 — )• 0, {X^,Y^) convergence to 
{X^,YI^) in the sense that 

vtG[o,r], E[|x7-x°|2 + |y,T-y,0|2] ^0. (6.2) 

The proof of Theorem 6.1 is similar to that of the Theorem 3.1; so, we omit it. 
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